Fidelity of fermionic atom-number states subjected to tunneling decay 
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Atom-number states are a valuable resource for ultracold chemistry, atom interferometry and 
quantum information processing. Recent experiments have achieved their deterministic preparation 
in trapped few-fermion systems. We analyze the tunneling decay of these states, both in terms of 
the survival probability, and the non-escape probability, which can be extracted from measurements 
of the full counting statistics. At short times, the survival probability exhibits deviations from the 
exponential law. The decay is governed by the multi-particle Zeno time which exhibits a signature 
of quantum statistics and contact interactions. The subsequent exponential regime governs most of 
the dynamics, and we provide accurate analytical expressions for the associated decay rates. Both 
dynamical regimes are illustrated in a realistic model. Finally, a global picture of multi-particle 
quantum decay is presented. 
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I. INTRODUCTION 

Successful production and preservation of atomic 
states containing an exactly known number of parti- 
cles (so-called atomic Fock or atom-number states) is 
of importance for ultracold chemistry, atom interferom- 
etry, quantum information processing and investigating 
the foundations of quantum physics. Among the meth- 
ods proposed for creating atom-number states are atom 
culling in time dependent traps [ - ] and the use of sim- 
ilar techniques in optical lattices [7, 8]. Spectacular ex- 
perimental progress in the deterministic preparation of 
few-atom states has been reported in [9, 10]. But when 
the confining potential has finite potential barriers, which 
allow tunneling leakage, an initial trapped state with a 
well-defined atom number eventually evolves into a mix- 
ture of several atom- number states, as shown in Fig. 1. 
This process can be quantitatively characterized by the 
fidelity decay. The fidelity is often defined as the proba- 
bility to persist in the initial atom-number state or, some- 
what less stringently, by the non-escape probability, i.e., 
the probability that the total number of atoms within 
the trap remains the same. Unlike the integrated density 
profile of an atomic cloud [11, 12], both probabilities refer 
to all particles occupying a certain subspace of a Hilbert 
space simultaneously and are, in this sense, multi-particle 
observables. 

The purpose of this paper is to investigate the fidelity 
decay of an atomic Fock state due to tunneling leak- 
age. Hence, it is both of fundamental interest to un- 
derstand the quantum decay of multi-particle systems, 
as well as of practical relevance to determine the de- 
cay rates of trapped Fock states [1-6, 9, 10]. We focus 
on ultracold atomic vapors confined in tight, effectively 
one-dimensional, waveguides and, more specifically, on 
polarized fermions and related fermionized systems such 
as a bosonic cloud in the Tonks-Girardeau regime [13]. 
Since the rate of atom losses due to three-body collisions 
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FIG. 1. Fidelity decay of an atom-number state of polarized 
fermions induced by tunneling losses. The probability to find 
n atoms inside the trap, p(n,t), is initially peaked at n = 
8. As the time of evolution goes by, p(n, t) broadens and 
shifts to lower values of the mean atom number. The trapping 
potential model is further discussed in Sect. VI. 



is suppressed due to spatial anti-bunching [14, 15], these 
systems are well suited for investigating the quantum dy- 
namics of multi-particle tunneling decay. Although the 
tunneling decay is expected to roughly follow an exponen- 
tial law, deviations are expected both at short and long 
times on theoretical grounds. Long-time deviations from 
exponential decay in the tunneling dynamics of multi- 
particle systems have been studied in [16], and we shall 
focus on the short-time and exponential regime. Addi- 
tional relevant works include studies of the fidelity decay 
in a multi-particle Loschmidt echo [17, 18] and different 
scenarios of two-particle quantum decay [19-23]. The 
rest of the paper is organized as follows: in Section II 
we define the non-escape and survival probabilities for a 
general quantum system. In Section III we analyze their 
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short-time evolution and the Zeno effect [24-20]. In Sec- 
tion IV we discuss the decay of one-dimensional trapped 
fermions and related systems and investigate the effect 
of quantum statistics in terms of the multi-particle Zeno 
time. Section V is devoted to the regime governed by an 
exponential decay law, an accurate semiclassical approx- 
imation for the corresponding decay rates is provided. 
Section VI illustrates the different dynamical regimes in 
a trap model of relevance to recent experiments [9, 10]. 
A general picture of multi-particle quantum decay is out- 
lined in Section VII, before closing the paper with a set 
of conclusions in Section VIII. 



II. NON-ESCAPE AND SURVIVAL 
PROBABILITIES 

We start by considering the time evolution of the prob- 
ability that a multi-particle quantum system, initially 
prepared in a (possibly mixed) state go, will be found 
inside a certain sub-space T> C H of the Hilbert space H. 
The ability to experimentally measure the Full Count- 
ing Statistics (FCS) [9] motivates the introduction of 
the multi-particle non-escape probability [16]. Consider 
the full atom-number distribution p(n,t) = (5(nx> — n)) t 
where n-p is the atom-number operator in the subspace 
of interest T>, n = 1, . . . , N is an integer, and the ex- 
pectation value is taken with respect to a time evolving 
state. p(n, t) represents the probability of finding exactly 
n-particles at time t in the subspace T>. Its typical be- 
havior is exhibited in Fig. 1 for a N-particle metastable 
Fock state. The probability to preserve the iV-particle 
Fock state at time t is given by p(n — N, t) and we shall 
denote it by Pv(t) in the following. 

Letting II be the projector on T>, we write this proba- 
bility as 

V v (t) = tr[e n(*)], ft(t) = exp(iHt)flexp(-iHt). (1) 

Let \n) be a basis on Ti, so that an initial state go — 
J2 n Pn\n)(n\. Let A = J2 n \ n )( n \ ^ e tne projector associ- 
ated with the space spanned by the initial state go, such 
that Ago A = go- The behavior of Px>{t) largely depends 
on whether go is contained within T>. If it is, II can be 
decomposed as the sum of projectors 

n = A + Q, Qg o Q = 0. (2) 

In this case we will refer to "Px>(i), satisfying 7 7 x>(0) = 1, 
as the non-escape probability probability P(i) which must 
be a non-increasing function at least for short times. In 
the special case of Q — 0, Pv(t) becomes the survival 
probability, i.e., that for the system to remain in its ini- 
tial state. This we will denote by S(t), bearing in mind 
that S(0) = 1. Note that P(i) > S(t) since the condition 
for the system to stay within V, is less strict than that 
for staying in one specific state go- An obvious exam- 
ple is a wavepacket contained within a spatial region A, 
for which the non-escape condition means not leaving A, 



as opposed to remaining in the initial state (survival). 
Finally, if go is not contained in T> we have 

go = flg U + Sg (3) 

where Sgo is the component of the initial state orthogo- 
nal to T>. In this case Pv(t) ma y increase even at short 
times, a simple example being a wavepacket arriving in 
an initially empty region of space. Next we consider the 
short time evolution of the two probabilities. 

III. THE SHORT TIME LIMIT 

The short time expansion of the probability to remain 
in the sub-space V takes the form 

P©(t)=tr( eo n)-ttr(ft)[n,fl])i 

-tr[g ( ^ {4 H 2 } - HUH)]t 2 + 0(t 3 ), 

(4) 

where {A,B} denotes the anticommutator of A and B. 
As it is well known (see, e.g., [24, 2- r )]), vanishing of the 
term linear in t leads to the Zeno effect: frequent checks 
of whether the system is still contained in T> (e.g., projec- 
tion measurements of II) would prevent it from leaving 
the sub-space. In particular, with the interval between 
checks being r, Vx>{t) would typically exhibit an effective 
exponential decay 

Vv(t) ~ exp(-Tt), i = tIt 2 z , (5) 

where the Zeno time tz, is determined by the coefficient 
multiplying t 2 in Eq. (4), 

rz = [tr[eo({n, H 2 }/2 - HUH)}- 1 ' 2 . (6) 

Nonetheless, we shall be concerned with quantum decay 
in the absence of measurements. It is readily seen that 
the linear in time term in Eq. (4) vanishes provided the 
initial state is contained in T>, i.e. the orthogonal compo- 
nent 5go — 0, or can be neglected, \\5qo\\ <§; 1. The decay 
of Vx>{t) becomes then quadratic in time and governed 
by tz according to 

Vv(t) = l-(t/rz) 2 + 0(t 3 ). (7) 

Under this condition, from Eq. (G) for the Zeno time, we 
have 

T Z = [AH go - tv(g Q HQH)]- 1/2 , (8) 

and not merely given by the inverse of the energy variance 
of the initial state, A_ff eo = tr[goH 2 ] — tr[goH] 2 . The last 
term in the square bracket vanishes in the case of the 
survival probability, Q — 0. Thus, a survival Zeno time 
never exceeds that in the non-escape case, reflecting the 
fact that it is easier to maintain a system within a larger 
subset of its Hilbert space. (Note that tz in Eq. (8) 
becomes infinite, as it should, if T> is chosen to coincide 
with H). Next we proceed to the case of several particles 
confined in a potential trap with tunneling leakage. 
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IV. MULTI-PARTICLE ZENO TIMES OF 
FERMIONIZED SYSTEMS 

Consider N non-interacting particles initially trapped 
in a potential well Vj(x) which, at t = 0, is instantly 
converted into a trapping potential with a finite barrier, 
V(x), as that in Fig. 2 below, thus allowing the particles 
to escape into the continuum. The Hamiltonian of the 
system reads 

N 

H(xi, ...,x N ) = ^ h{xi), 

h( Xi ) = -d 2 x j2 + V(xi), i = 1,2, N. (9) 

Spin polarized fermions fall within this description, and 
certain strongly interacting systems can be described in 
a similar way. This is the case of bosonic atoms in 
the Tonks- Girardeau regime, where strong zero-range 
hard-core interaction leads to fermionization [ ]. In 
this case there is a one-to-one correspondence (the 
Bose- Fermi mapping [13]) between the symmetric state 
of strongly interacting bosons, \l/ TG (xi, . . . , in), and 
the anti-symmetric state of the dual system of non- 
interacting fermions, ^ F (xi, . . . , xn), 

*tg(xi, . . . ,x N ) = A^> F (xi, ■ ■ -,x N ), (10) 

where the antisymmetric unit function A = 
J| 1<J<fc<N sgn(xfc — Xj). It follows that for the 
Tonks- Girardeau gas both the non-escape and the 
survival probabilities coincide with those calculated for 
the corresponding system of non-interacting fermions 
[16], 

P TG( t ) = pf ( t)) S TG (<) = S F (i). (11) 

Both of these systems were found to be optimal for the 
preparation of atomic Fock States using atom culling 
techniques [1-6, 9, 10]. In the following discussion we 
will simply refer to a "fermionic" system, bearing in mind 
that the results apply both to the Tonks- Girardeau gas 
with infinitely strong contact interactions and to a sys- 
tem of noninteracting polarized fermions. 

The ground state of a fermionic system is given by the 
Slater determinant, and we choose 

<(xi,...,x N ) = -^det* fe=1 [0„(x fe )], (12) 

where \4> n ), n = 1, •••jN, are the N lowest eigenstates of 
the one-particle Hamiltonian, h\<j) n ) = £ n \4>n)- In order 
to study the effect of quantum statistics on the Zeno 
time it is instructive to consider also the ground state for 
distinguishable particles, 

N 

^ is '(xi,...,X N ) = Y[<p n (x n ), (13) 
n=l 



and the somewhat artificial case of a bosonic excited 
state where the same first N levels are occupied by non- 
interacting particles obeying Bose-Einstein statistics, 

*f (xi,...,x N ) = -^=per^ fe=1 [0 n (x fe )], (14) 

where per stands for the permanent, i.e., the sum in the 
r.h.s. of all permutations of indices k for a fixed order of 
indices n. Note that these three states, Eqs. (12), (13), 
and (14), have all the same energy. 

The survival Zeno times [Eq. (8) with Q = 0] for the 
three cases are calculated by inserting in Eq. (8) the 
appropriate initial state, (12), (13) or (14). For distin- 
guishable particles the calculation is straightforward, for 
polarized fermions the matrix elements in Eq. (8) can 
be evaluated using the Slater-Condon rules [27], and for 
bosons an extension of the later is required to symmetric 
states. The result can be written in a compact form, 

rz = lj2 Ak - + 2a T,\^\ h \^\ 2 } 2 ' ( 15 ) 

^ n n <k ' 

where Ah n = (<p n \h 2 \<p n ) - (4>n\h\<l>n) 2 , and a = for 
distinguishable particles, a = — 1 for fermions, and a = 1 
for excited bosons. 

In Eq. (15), the first sum, common to all statistics, 
depends on the spread of the energy in the initial one- 
particle states induced by the tunneling decay. This con- 
tribution dominates the Zeno time. The second sum, 
arising from the indistinguishability of the particles and 
the associated symmetrization of the initial state, and 
more precisely, from its immanant (either determinant or 
permanent) structure. While it reduces the survival Zeno 
time for non-interacting excited bosons, it increases Tz 
in polarized fermions and bosons in the TG regime, lead- 
ing to a slowing down of their decay. The symmetriza- 
tion imposed by exchange quantum statistics plays no 
role, as it follows from the fact that the multi-particle 
Zeno time is shared by dual systems related by the Bose- 
Fermi mapping. The different corrections to Tz arise only 
for indistinguishable particles, and are manisfest thanks 
to the effect that contact interactions (including as such 
the Pauli exclusion principle) have on the energy disper- 
sion of the initial state. We recall that short time decay 
V v {t) = 1 - {t/rz) 2 = 0(t 3 ) is governed by the Zeno 
time. Also, we recall that when the system is frequently 
observed its lifetime 7 _1 becomes proportional to r§, see 
Eq. (5). It follows that a fermionic state (or bosons in 
the TG regime) decays at short-times more slowly than 
non-interacting excited bosons. This result is somewhat 
counter intuitive in nature, given that as a result of the 
Pauli exclusion principle (or hard-core contact interac- 
tions in the TG gas) fermionized systems exhibit spa- 
tial anti-bunching, while non-interacting bosons prefer 
to group together (bunching). However, this intuition 
applies only in subsequent stages of evolution, while the 
short-time dynamics is exclusively governed by the en- 
ergy dispersion of the initial state. 
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In the following section we shall describe the rates asso- 
ciated with the subsequent exponential decay, where the 
indistinguishability of the particles and density-density 
correlations play a subdominant role with respect to 
the energy distribution. The intuition based on spatial 
bunching applies to the long-time asymptotics of multi- 
particle quantum decay. At long times, deviations from 
the exponential decay occur due to the possibility that 
the decay products recombine to reconstruct the initial 
state. As a result, spatial bunching and antibunching 
effects play the dominant role, and fcrmionized systems 
decay according to a power law VT>(t) oc l/t a with an 
exponent a > about N times larger than in the case of 
non- interacting bosons [ ] . 
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FIG. 2. The 'bathtub' trapping potential V(x) (solid line). 
Also shown is the initial well Vi(x) (dashed line) 



V. EXPONENTIAL REGIME AND DECAY 
RATES OF FERMIONIZED SYSTEMS 

We shall now turn our attention to the regime char- 
acterized by exponential decay, most easily observed in 
experiments. While manipulating cold atoms one may be 
interested in maintaining exactly N of them in a specified 
region of space A, usually large enough as to enclose the 
trap subspace. The corresponding non-escape probabil- 
ity is obtained from Eq. (1) by choosing the projector 
onto the corresponding part of the multi-configurational 
space, 



ft = Y]_ / \Xn)(Xn\dx n . 



(16) 



For a fermionic system starting in its ground state (12) 
the non-escape probability reads [16] (we drop the super- 
script F in the following and include a subscript for the 
particle number N) 



P N (t)=det£ fc=1 [^„(t)| XA |<M*)^ 



(17) 



where ((f> n (t)\xA\(f>k(t)} = J A dx4>* n {x, t)4> k {x, t), while for 
the probability to survive in the state (12), S^(t), one has 
[16-18] 



S N (i) = detl k=1 [(M0)\Mt))} 



(18) 



Once the initial state is prepared, after the Zeno time 
scale, each of the single-particle states {!</>&)} involved 
in the Slater determinant (12) experiences a nearly ex- 
ponential decay, i.e. \<p k (t)) » exp(-7 fc t)||</> fc (0)), k = 
1, . . . ,N. Then, the N-particle survival and non-escape 
probabilities decay exponentially, 



P N (t) 
Sn(*) 

with a decay rate 



cxp(-r<)P N (0), 

cxp(-r<)s N (o), 



r = 2£ 



(19) 



(20) 



Individual decay constants "f k are given by the imaginary 
parts of the complex energies E k corresponding to the 
first N leading poles, q k , of the reflection amplitude for 
the potential trap in the complex momentum q— plane, 
7fe = lmE(q k ) — Keq k lmq k . One can avoid precise de- 
termination of the momentum pole positions by making 
instead a simple semiclassical estimate [28]. Let E k be 
the energy of the fc-th bound state in the initial trap, and 
Xq < x\ < x\ the three turning points in the quenched 
potential satisfying V(x^) — E k (i = 0,1,2), e.g. for 
a potential where the left barrier is arbitrarily wide so 
that tunneling losses occur only through the right bar- 
rier, as in the case depicted in Fig. 2. The semiclassical 
probability to tunnel in one attempt across the potential 
barrier is given by (see e.g, [29]) 



T(E k ) 



exp[-25( J B fe )] 



{l+cxp[-2S(£,)]/4} 2 - 



(21) 



where S(E k ) = Q &x{2^/(x) - E k \Y 12 is the complex 
action corresponding to the classically forbidden region 
x\ < x < X2- Since the period of the bound motion is 

given by r k = 2 dx{2[V{x) - £ fc ]}- 1/2 , the particle 
impacts on the barrier with an approximate frequency 
n k = 1/ r k . Multiplying T by the number of impacts per 
unit time gives individual decay rates 



27fc 



l T{E k ), 



1, 



,N, 



(22) 



which together with Eq. (20) yield the desired rate of 
the exponential decay in Eqs. (19). 



VI. THE MODEL AND RESULTS 

As a realistic model of a one-dimensional trap we con- 
sider a smooth bathtub potential (Fig. 2, dashed) 



fc=i 



Vj{x) = ~V 



1 — tanh 



\x\ -L/2 



(23) 
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FIG. 3. Non-escape and survival probabilities for fermionic 
atom-number states with N = 4 and N = 8 leaking out of 
potential (24) with C — 8. The insets show the quadratic 
behavior of Pn(£) and Sjsr(i) in the short-time regime (solid) 
and the fitted quadratic parabola (dots). 



whose right wall is instantly turned into a barrier of finite 
width at t = (Fig. 2, solid), 



V{x,t) = ~V 



1 — tanh 



- L/2 



Q(a- x) 



1 + tanh 



x — a — Li/2 

0\ 



&(x- a). 



(24) 



Here Vb is the depth of the initial well (and also the bar- 
rier height of the metastable potential), L and L\ are the 
widths of the well and barrier, respectively, and a ipx) 
determine the smoothness of the inner (outer) potential 
walls. It is convenient to introduce the dimensionless 
variables (we re-introduce the Plank's constant K) 

x^x/L, t-tt/to, V ->V t /h (25) 
with to = mL jh. The absorbing potential introduced 
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FIG. 4. Non-escape (a) and survival (b) probabilities for a 
fermionic 12-atom state leaking out of potential trap in Eq. 
(24) with capacity C = 12 (solid) (all times are in units of 
to). Also shown are the short time quadratic fits (dashed) 
and the times t q beyond which the quadratic approximation 



fails (vertical dashed). 



by Manolopoulos [30, 31] is employed to avoid unphys- 
ical reflections at the boundaries of the numerical grid, 
L\°2 ■ Bound one-particle eigenstates of the initial well, 
<p n (x,0), are obtained by a standard finite-difference 
technique, and then evolved in time using the Crank- 
Nicholson scheme, to yield <j> n (x,i) required in Eqs. (17) 
and (18). We use L x /L = 0.08, a/L = a x /L = 0.01, 
a/L = 0.55 and V t /h = C 2 n 2 , where C is the capacity 
of the well, i.e., the maximum number of bound states 
it supports. We also note that for 23 Na atoms in a po- 
tential well with L = 80 ^m to is about 2.39 s. Finally, 
we chose L b ° x / L = -20, L\ ox / L = 30, and the absorbing 
potential identical to that used in [G] . For the calculation 
of the non-escape probability (1) for a fermionized sys- 
tem in its ground state (12) we chose the spatial region 
A in Eq. (16) to include most of the trap, 



A = (-oo, a], 



(26) 



so that Sqq in Eq. (3), now associated with exponential 
tails of the (j> n {x) extending into the right classically for- 
bidden region of Vi(x), is small. Figure 3 shows the de- 
cay dynamics of Pn and Sn for different fermionic atom- 
number states when the initial well supports a maximum 
of C = 8 bound states. The two curves are remarkably 
similar, given that Sn is sensitive to the population of 
individual one-particle states, while Pn only depends on 
whether atoms have left the region A. In agreement with 
the discussion in the previous sections, two regimes can 
be identified. The short-time behavior is characterized 
by a decay quadratic in time as shown in the inset of 
Fig. 3. Subsequently, as dictated by Eq. (19), the expo- 
nential law holds. The associated transition is illustrated 
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FIG. 5. Dependence of the survival Zeno time tz on the num- 
ber of trapped fermions N for a potential (24) with (7 = 8. 
The analytical Zeno times given by Eq. (15) (squares) are 
shown together with the numerical values extracted from a 
parabolic fit to the exact short-time decay dynamics (dots). 
The inset shows the corresponding exponential fits to the de- 
pendence of tz on N. 



in Fig. 4 for a N = C = 12 fermionic Fock state, and 
occurs approximately in the time scale 



r 2 F 



(27) 



where the Zeno time tz and decay rate T are given by 
equations (15) and (20), respectively. 

We first focus on the short-time dynamics, where the 
characteristic scale of the decay is given by the Zeno time 
tz- As the energy dispersion of the initial state with re- 
spect to the Hamiltonian for t > increases, the corre- 
sponding tz is reduced. In particular, this is expected 
from Eq. 15 for increasing particle number N. Figure 
5 shows the survival Zeno time tz as a function of N. 
Its dependence on the particle number is enhanced as N 
approaches C and the initial Fock state involves energy 
components closer to the brim of the trap, which result in 
the spatial extension of the atomic cloud beyond the trap 
region A] due to existence of tunneling tails for x > a. 
Indeed, it is found numerically that the suppression of 
the Zeno-time is approximately exponential with respect 
to N, 7-2 (N) « 7-2(1) exp(-N/2), see the inset in Fig. 5. 

This rapid decline of rz(N) might suggest that the 
short-time deviations from the exponential law should 
be irrelevant for Fock states with already moderate num- 
bers of atoms. However, this is not completely the case, 
as is illustrated in Fig. 4 which demonstrates that the 
fidelity decay of an Fock state also accelerates with the 
number of atoms, and becomes significant already at the 
end of the quadratic evolution, provided N is sufficiently 
close to C. To quantify this effect, we evaluated the tran- 
sition time t g (N) (Fig. 4, vertical dashed) such that for 
t > t q the quadratic approximation to Pn(£) or Sn(*) 
starts to fail. (Note that t q is not the same as tz, which 
determines the curvature of the short-time parabolic de- 
cay.) The plot of the fidelity factor P N (i g )/P N (0) vs. N 
shown in Fig. 6 confirms that as N/C approaches unity 
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FIG. 6. Loss of fidelity for fermionic atom- number states 
leaking out of potential trap described by Eq. (24) for two 
different capacities C = 12,30 (N < C). The fidelity is eval- 
uated at the transition time t q beyond which the quadratic 
approximation fails (cf. Fig. 4), this is, right before the ex- 
ponential regime sets in. 
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FIG. 7. Dependence of the exponential decay rates for Sn(£) 
(square) and Pn(£) (triangle) on the number of trapped atoms 
N for a potential trap in Eq. (24) of capacity C = 8. The 
observed survival decay rates are in excellent agrement with 
the analytical prediction based on Eqs. (22). 



a significant fidelity decay of an atom number state oc- 
curs before the exponential regime sets in. The effect 
is more pronounced for small traps, the goal of current 
experimental efforts. 

We next focus on the characterization of the exponen- 
tial regime. The dependence of the decay rate (20) on 
the number of particles N is in good agreement with the 
semiclassical estimate derived from Eqs. (22), as shown 
in Fig. 7. This is to be contrasted with the case of a trap 
with thin potential barriers [11] where the decay rate of 
any single-particle excited state is soon governed by the 
longest-lived resonance, favoring the observation of mul- 
tiple exponential regimes (for N = 2 with a <5-barrier, see 
[21]). In such case the multiparticle decay rate becomes 
r = 2jiN linear in N. Note that the expression Eqs. (20) 
for the decay rates equally applies to distinguishable par- 
ticles or excited states of non-interacting bosons. Its suc- 
cess in reproducing the decay rates extracted from expo- 
nential fits to the exact decay dynamics, points out that 
the indistinguishability of the particles and the many- 
body correlations among them play a subdominant role 
in the exponential regime. Nonetheless, as the initial 
state includes states closer to the brim of the trap, Eq. 
(21) losses accuracy. 
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VII. DISCUSSION 

Given the recent outburst of works dealing with few- 
body quantum decay [11, 12, 16-23], it is interesting 
to establish the current understanding and point out 
some aspects deserving further studies. In particular, 
combining the results obtained in this manuscript with 
those in [16], we can identified the following stages in 
multi-particle quantum decay of systems with contact- 
interactions. As in the single-particle decay, three 
regimes are found: 

• Short time asymptotics: Zeno regime. - A 
quadratic-in-time quantum decay occurs that is 
characterized by the multi-particle Zeno time, i.e. 
the inverse of the energy variance of the quenched 
Hamiltonian (that, for t > 0) evaluated in the ini- 
tial state, Eq. (15). The main contribution to 
this time scale arises from the energy spread of 
the initial state. Nonetheless, in the presence of 
interactions there is a counterintuitive correction 
reflecting the indistinguishability of the particles. 
This correction slows down the decay of polarized 
fermions with respect to degenerated states of non- 
interacting bosons. It would be clarifying to con- 
sider systems with Generalized Exclusion Statistics 
[32] and compare the multi-particle quantum decay 
between systems with finite-interactions related by 
Bose- Fermi duality [33]. 

• Exponential regime. - Following the Zeno regime, 
exponential decay sets in, characterized by multi- 
particle decay rates well described in a semi- 
classical approximation, i.e. Eq. (22). The de- 
pendence of the decay rates is dominated by the 
energy spread of the initial state. This is expected 
to be the general behavior for smooth potentials, 
but for finite-interactions the question remains as 
to the validity of the semi-classical approximation. 
For example, considering a ID Bose gas [22], does 
Eq. (22) evaluated at the quasi-momenta of the 
initial state (Bethe roots) describe accurately the 
decay rate, or are there important corrections aris- 
ing from correlations among the particles? Inter- 
estingly enough, the exact decay dynamics of two 
fermions for arbitrarily thin potential barriers ex- 
hibits a transition among different decay rates [21]. 
The observability of multiple exponential regimes 
for arbitrary particle number N and as a function 
of the potential features well deserves further stud- 
ies. 

• Long-time asymptotics: post-exponential power-law 
decay. - The subsequent regime is characterized 
by a power-law decay and arises due to the possi- 
bility that the decay products recombine to form 
the initial state in a classical sense. This regime 
is governed by the short-range density-density cor- 
relations, that lead to a dramatic dependence of 



the power-law exponents on the particle-number 
N [16]. The current understanding would benefit 
from studies for systems with finite and long-range 
interparticle interactions as well as two and three 
dimensional systems with degeneracy arising from 
angular momentum. The power-law exponents in 
these systems are expected to exhibit a rich depen- 
dence on N. 

We note that a deep understanding of multi-particle 
quantum decay is relevant across very different fields 
of physics, beyond quantum foundations. For instance, 
a good control of the particle-number might be key to 
advance the field of quantum simulation with ultracold 
gases and this goal is being currently pursued in differ- 
ent laboratories [9, 10]. Similarly, the short-time quan- 
tum decay is at the heart of certain applications such as 
dynamical decoupling schemes for high-fidelity quantum 
memories, quantum gates, quantum computing and gen- 
erally to extend the coherence times in quantum systems 
[ ]. And at the same time, the existence of different 
stages of quantum decay may have important cosmolog- 
ical implications [35]. 



VIII. CONCLUSIONS 

In summary, we have analyzed the many-body tun- 
neling decay of trapped fermionic atom-number (Fock) 
states in a realistic model of relevance to recent experi- 
ments [9, 10] . We focused our attention on the fidelity de- 
cay resulting from tunneling leakage, quantified as a de- 
crease in the non-escape or survival probability, and iden- 
tified the signatures of contact interactions and quantum 
statistics in the short-time multi-particle decay. Even 
though the survival and non-escape Zeno times rapidly 
decrease with the number of atoms N, as N approaches 
the maximum capacity of the trap a substantial loss of fi- 
delity is found already at the end of the quadratic (Zeno) 
evolution, before the exponential regime sets in. More- 
over, explicit expressions for the decay rates in the ex- 
ponential regime have been provided. Our results are 
amenable to experimental verification by standard tech- 
niques, e.g., by registering the time-evolution of the full- 
counting statistics in a trap which allows leaking by quan- 
tum tunneling (see, for example, Ref. [9]). Finally, we 
note that strong s-wave scattering in ultracold bosons or 
spin-polarization in fermions suppress three body-losses, 
facilitating the study of genuinely quantum losses in these 
systems. 
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